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We investigate theoretically the resonant inelastic x-ray scattering (RIXS) at the K edge of 
Si on the basis of an ab initio calculation. We calculate the RIXS spectra with systematically 
varying transfered-momenta, incident-photon energy and incident-photon polarization. We 
confirm the anisotropy of the experimental spectra by Y. Ma et al. (Phys. Rev. Lett. 74, 478 
(1995)), providing a quantitative explanation of the spectra. 
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1. Introduction 

Inelastic x-ray scattering is a promising method to study electronic structures in matters. 
It is advantageous to use a resonant enhancement by tuning photon energy near absorption 
edge. Recently it has been revealed that the resonant inelastic x-ray scattering (RIXS) is 
a powerful tool for elucidating the electronic properties of solids. RIXS measurement has 
been applied to search for charge excitations in several cuprates 1,2 and manganites, 3 with 
tuning incident photon energy to the absorption K edge. Since corresponding photon energies 
are in the hard-X-ray region, the momentum transfer can not be neglected. Analyzing the 
transfered-momentum dependence of the RIXS spectra, we can clarify the characteristics of 
charge excitations in strongly-correlated electron systems. 

In semiconductors such as Si and Ge, the interpretation of the RIXS spectra may be 
simpler because of weak electron correlations. The final state of RIXS consists of one electron 
in the conduction band and one hole in the valence band. With neglecting the interaction 
between the electron and the hole, one can draw useful information of single-particle spectra 
from the transferred-momentum dependence. Actually Ma has tried to determine valence band 
structures of Si by systematically varying scattering vectors in the RIXS experiment. 

We have already carried out a RIXS experiment at the K edge of Ge with the same aim. 
Although we systematically varied transfered-momenta, we obtained merely a broad inelastic 
peak as a function of photon energy, whose shape was nearly independent of transfered- 

* E-mail address:nisikawa@spring8. or.jp 



1/?? 



J. Phys. Soc. Jpn. Full Paper 

momenta. Performing the band structure calculation, we obtain the spectral shape in agree- 
ment with the experiment, and clarified the origin of the spectral shape. 5 This analysis led 
us to realize that the decisive factor whether the spectra are depending on the transferred 
momenta or not is the core-level width. If the core-hole level in the intermediate state is sharp 
enough, we can definitely select the momentum of the excited electron in the conduction band 
by sharply tuning the incident-photon energy, and thereby we can also specify the momentum 
of the hole in the valence band by setting the scattering geometry. In this situation, sharp 
peaks are expected to be observed as a function of scattered photon energy. If the core-level 
width is not small, many channels of exciting an electron with different momenta in the con- 
duction band are expected to be opened. Thereby, sharp peaks are blurred up by overlapping 
each contribution. In an extreme situation, the spectra become close to a superposition of 
the density of states (DOS) projected onto p-symmetric states. Unfortunately, for Ge, the Is 
core-level width is as large as ~ 2eV, an order of the width of conduction band, so that the 
RIXS spectra did not show clear momentum dependence. 

For Si, since the K edge energy is around 1840 eV, the transfered-momentum cannot 
be neglected in the RIXS experiment. The Is core-level width is estimated about 0.6 eV, 
which is smaller than the width of conduction band of Si. These facts imply a possibility of 
observing the transfered-momentum dependence at the K edge in the RIXS spectra. Actually, 
Ma et. al observed such transfered-momentum dependence. 6 In this paper, we analyze Ma's 
experimental RIXS spectra with the help of ab initio band structure calculation, and elucidate 
the relation between the spectral shape and the underlying electronic structure. 

This paper is organized as follows. In §2, the formalism for calculating the spectra is 
described. In §3, the calculated results and discussion are presented. Section 4 is devoted to 
concluding remarks. 

2. Formulation 

RIXS is described by a second-order process that the incident photon with energy Hloi, 
momentum ftx\ 1 and polarization ei is virtually absorbed by exciting the core electron to the 
conduction band and then a photon with energy Huj2, momentum /iq2 and polarization e2 
is emitted by filling the core-hole state with a valence electron. The scattering geometry is 
shown in Fig. 1. The normal axis n and theta axis t in this figure are respectively defined as 
follows: n oc q2/|q2| — qi/|qi|, t oc qi x q2. The setup for scattering is uniquely determined 
by choosing the normal and theta axis from crystal axes of Si. 

To obtain the double differential scattering cross-section, we use the generalized Fermi's 
golden rule where the interaction between photon and electrons is treated by second order 
perturbation theory. An independent-particle approximation seems to be appropriate, since 
electron correlations are expected to be weak. 4 In such a situation, the double differential 
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where /iw = /icji — k<j2, and Till = Tiqi — 7iq2 are energy and momentum of the final state. 
e e (k), e^(k') and e c are the energy of the excited electron with crystal momentum k in the 
conduction band e, that of the hole with crystal momentum k' in the valence band h, and that 
of the core state, respectively. Overlined quantity indicates the complex conjugate. The F is 
the Is core-level width of Si. Quantity t a (k',/i|e2) = / drVv,fe( r ) e 2 ■ P</>i s (r — R-a) describes 
the transition from the valence band to the Is core state, where R a , and </>* s are the 

position vector of atom a in unit cell, Bloch-wave function of an electron in the valence band 
h with crystal momentum k', and ls-atomic orbital, respectively. The crystal momentum 
conservation for the whole process is contained in the factor of Kronecker S, 



where G is the set of reciprocal lattice vectors. 
3. Calculated Results and Discussion 

3. 1 Band Structure Calculation of Si 

We perform a band structure calculation using the full-potential linearized augmented- 
plane-wave (FLAPW) method within the local-density approximation (LDA). The local 
exchange-correlation functional of Vosko, Wilk and Nusair is employed. The angular mo- 
mentum in the spherical-wave expansion is truncated at £ max = 6 and 7 for the potential 
and wave function, respectively. The energy cutoff of the plane wave is 12 Ry for the wave 
function. Figure 2 shows the energy vs. momentum relation thus evaluated. The energy band 
is labeled by attached numbers for later use. 

3.2 RIXS spectra 

3.2.1 Anisotropy of RIXS at K edge of Si 

We first analyze the transfered-momentum dependence of RIXS spectra for the incident 
photon energy TujJ\ — £min — where e m i n is the minimum energy of conduction-bands. When 
T is very small, the denominator in Eq. 1 forces the excited electron stay just at the bottom 
of the conduction band. There are six such points, which are denoted as kj (i = 1, ■ • ■ , 6) in 
Fig. 3. They belong to band 5 (see Fig.2). We can specify the momentum of hole = k« — fl in 
the valence band by specifying the scattering vector. In such a situation, the double differential 
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scattering cross-section can be rewritten as 
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The sharp peaks are expected at Huj = hio(i,h) = e m j n — e^(qj) in the spectra, where h 
(= 1, • • • ,4) indicates the valence band index. When the scattering vector is along a sym- 
metry line, the number of peaks is reduced because of degeneracy of peaks and no inten- 
sities to some peaks. For n = (100), for example, the momenta of the holes are shown in 
Fig.3, and the corresponding valence band energies are the same due to symmetry, that 
is, e/j(q 2 ) = e h (q 3 ) = e fe (q 4 ) = e h (q 5 ) with h = l,---4, except for qi and q 6 . Therefore, 
several excitation energies coincide with each other, and thereby different peak positions 
are huj = hu}(l,h),hw(2,h),hu}(6,h) with h = 1, •■■ ,4. In addition, the spectral weights 
A(i, h\£l, ei = a) for momentum ki and kg are found to be almost zero. Thus only four peaks 
appear at Hu = hu>(2, h)(h = 1, ■ • • , 4), as shown in Fig. 4 (c). With slight increase of the T 
value, the excited electron can stay at several states near the bottom of the conduction band, 
and thereby the valence hole can stay at several states near qj(z = 2, • • • ,5). Therefore the 
spectral peaks are obscured. We can make similar analysis for n = (111), (110). Such spectral 
changes are demonstrated in Fig.4 for n = (111), (110) and (100), with F changed from 0.02 
eV to 0.6 eV. 

When T is large, the momentum of the excited electron can take whole values in the first 
Brillouin zone, and so does the momentum of the hole in the valence band. The spectral 
peaks are smeared out, and thereby becoming close to the DOS projected onto p symmetric 
states. Thus we have the spectra almost independent of scattering vector. This behavior is 
demonstrated in Fig. 5, where T = 1.6 eV. 

For the actual Si case, T is close to 0.6 eV, and the spectral shape is expected to show 
the transferred-momentum dependence. Figure 6 shows the RIXS spectra with T = 0.6 eV 
for huji = 1840.8 eV, in comparison with experimental results by Ma et. al. 6 In the experi- 
mental setup of Ma et. al 6 the normal axis n is (111), (110), and (100), but the theta axis 
t and incident-photon polarization are not clearly indicated in their paper. 6 Assuming the a 
polarization for the incident photon and the theta axis t as given in the figure, we are suc- 
cessful in reproducing the experimental shape such as three peak structure for n = (111) and 
the different spectral shape for n = (110) and (100). Corresponding to the above discussion, 
we can identify where the peak structures come from by decomposing the spectra into each 
contribution from the hole band and the electron band. Figure 7 shows such decomposition, 
where attached numbers indicates the band index given in Fig. 2. 
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3.2.2 Incident-photon energy dependence of RIXS spectra 

The spectra strongly depend on the incident-photon energy. Figure 8 shows the spectra 
with varying incident-photon energy from 1840.0 eV to 1841.2 eV. The scattering vector is 
fixed to be proportional to n = (111). The characteristic three peaks A, B and C are visible 
only in the narrow energy region between 1840.2 eV and 1840.8 eV. For Juoi = 1842.5 eV, the 
anisotropy almost disappears, as shown in Fig. 9. The result is consistent with experimental 
results of Ma et. al. 6 

3.2.3 Incident-photon polarization dependence 

Figure 10 shows the RIXS spectra calculated for a and tt incident polarizations. It is 
found that the difference in RIXS spectra in higher Huj2 for n = (100) is large, while not 
for n = (111) and (110). This tendency is more clearly seen in the limit situation T = 0.02 
eV where ftw\ = e m - m — e c , as shown in Fig. 11. The positions of peaks expected to be made 
by each electron- hole pairs are also presented in Fig. 11. To make the discussion clear, we 
consider the limit situation mentioned above. Then, the Is electron is excited to the states 
with momentum kj (i = 1, • • • , 6). Corresponding wave functions projected onto p-symmetric 
states centering on R a (a = 1, 2) are given by Pj=i^k lj5 oc p x , P; = i^k 2 ,5 oc p y , P;=i^k3,5 oc p y , 
P*=i^k 4 ,5 oc p z , P/=i^k 5 ,5 oc p z , Pi=iipk 6>5 oc p x , as shown in the left side of Fig.ll. In the 
case of n = (100), t = (011), the incident-photon polarization vector ei is proportional 
to (0,1,1) in the a polarization and (\/2, — 1, 1) in the tt polarization. The x-component 
of ei is not zero only in the case of tt polarization. Then the t a (kj,5|ei) in Eq.4 shows a 
large change at ki and kg, with changing polarization; \t a (ki : Q, 5\tt)\ » |i a (ki,6, 5|<r)| ~ 0, 
Ko(k2,-,5,5|7r)| ~ |t a (k2 ) ... ,5, 5|er)|. Therefore, we can expect that the electron-hole pairs with 
momenta ki, qi and k@, q6 contribute to RIXS spectrum in the case of the tt polarization. 
Note that the positions of the peaks made by the electron-hole pairs with momenta ki, qi 
and k 6 , q 6 in the tt polarization are quite different from these of the peaks made by other 
electron-hole pairs with momenta kj, qj(i = 2, • • • ,5). From Fig. 11(c), we find that electron- 
hole pairs with momenta ki, qi for h = 2, 3, 4, and kg, q6 for h = 3, 4 significantly contribute 
to forming the peak around IvjJi = 1839 eV in the tt polarization. Such a contribution gives 
rise to the appreciable intensity around Hu>2 = 1839 eV in the realistic RIXS spectrum for the 
tt polarization, as shown in Fig. 10. The same discussion mentioned above can be applied in 
the case of n = (111), t =(1-10) and n = (110), t =(1-10). In the case of n = (111), t =(1- 
10) and n = (110), t =(1-10), the incident-photon polarization vector ei is proportional to 
(1,-1,0) in the a polarization and (1,1, a), (a / 0) in the tt polarization. The z-component 
of ei is not zero only in the case of tt polarization. Therefore, the t a (kj,5|ei) shows a large 
change at k4 and ks, with changing the polarization; |i a (k4,5, 5|vr) | >> |i a 0*4,5 > 5|a)| ~ 0, 
|io(ki,2,3,6j 5 1 7r) I ~ |t a (ki,2,3,6j 5|ct)|. Therefore, we can expect that the electron-hole pairs with 
momenta k4, q4 and qs contribute to RIXS spectrum in the case of the tt polarization. 
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From Fig. 11(a), (b), we can see the following facts. In the case of n = (111), t =(1-10), the 
positions of the peaks made by the electron- hole pairs with momenta k4, q4 and ks, qs in the 
7r polarization are same as these of the peaks made by other electron-hole pairs with momenta 
kj, c[i(i = 1,2,3,6). In the case of n = (110), t =(1-10), the position of the peak made by 
the electron- hole pairs with momenta k4, q4 and ks, qs for h = 2 in the it polarization is 
same as that of the main peak in the a polarization. On the other hand, the positions of the 
peaks made by the electron-hole pairs with momenta k4, q4 and ks, qs for h = 3,4 in the 
7r polarization are different from that of the peaks made by other electron-hole pairs with 
momenta kj, qj(z = 1, 2, 3, 6). But the new peaks made by electron-hole pairs with momenta 
k4, q4 and ks, qs for h = 3,4 are not so significant. Summarizing, in the case of n = (111), 
t =(1-10) and n = (110), t =(1-10), significant peaks at new positions are not expected by 
varying incident-photon polarization. Therefore, the realistic RIXS spectra for n = (111), 
t =(1-10) and n = (110), t =(1-10) are not sensitive to the incident-photon polarization, as 
shown in Fig. 10. 

4. Concluding Remarks 

We have calculated the RIXS spectra at the K edge of Si on the basis of the band 
structure calculation. We have successfully reproduced the anisotropy of the RIXS spectra 
in the experiment by Ma et. a/. 6 We have analyzed the spectra with systematically varying 
transfered-momenta, incident-photon energy and incident-photon polarization, thus having 
provided a quantitative explanation of the spectra. With the help of the band calculation, one 
can draw information on the valence band structure from the RIXS spectra, although it is 
indirect. This information seems helpful with considering that the RIXS is bulk-sensitive in 
contrast to photoemission with a lower energy photon which is surface-sensitive. 

Finally, we comment the following general tendencies. To have an anisotropy in RIXS at 
K edge of specific element in semiconductor or band-insulator, the sizable momentum transfer 
and core-level width smaller than conduction-bandwidth are required. The photon energy is 
required to obtain sizable momentum transfer but corresponding core-level widths are usually 
large. Therefore, it is difficult to satisfy the requirement of the sizable momentum transfer 
and small core-level width. 
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Fig. 1. Geometrical relations between qi, q 2 and crystal axis of Si. 
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Fig. 3. Cross-section of the first Brillouin zone. The crystal momenta of electrons and holes generated 
in final state are presented. 
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Fig. 4. PJXS spectra for realistic value of core-level width (r = 0.6 eV) in comparison with the RIXS 
spectra for sufficiently small value of core-level width (r = 0.02 cV). 
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Fig. 5. PJXS spectra for large value of core- level width (r = 1.6 eV). The anisotropy almost disap- 
pears in the spectra. 
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Fig. 6. PJXS spectra obtained by (a)experiment 6 and (b)theory 
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Fig. 7. The spectra for hu\ =1840.8 eV, n=(lll), (110), and (100) are decomposed into each contri- 
bution of band-to-band transitions specified by band indices. 
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Fig. 8. Calculated FJXS spectra for n=(lll) with varying incident-photon energy between 1840.0 
eV and 1841.2 eV. 
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Fig. 9. Calculated RLXS spectra for tvjJi = 1842.5 eV with varying n=(lll), (110), and (100). 
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Fig. 10. Calculated PJXS spectra with varying incident-photon polarization. 
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Fig. 11. Incident-photon polarization dependence of RIXS spectra for hx>\ = e m i n — e c with sufficiently 
small value of r(=0.02 eV). The positions of peaks expected to be made by each electron-hole 
pairs are also presented. 
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